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1 Calculation of cell fate 

After infection of target cells, infected cells I\ can either become defectively infected (D), latently 
infected (Li), persistently infected (Mi) or activated virus-producing cells ih) or die during the 
intracellular eclipse phase. The proportions of cells that end up in a particular subpopulation 
were calculated as follows: 

n D-T- X T (1) 

8 h Ti 8 h r 2 

ttl=- x X- , (2) 

8 h Ti 8 h j 2 8 h j3 

Xm=t x — z x — z x- , (3) 

ShTi 5 hT2 S h r 3 8 h Y4 

Ti (1-/d)72 d-/i)73 (1-/m)74 7s 
7T/ = x x x x , (4) 

07x71 8 l2 j2 8 h j3 8 l4 j4 8 k y 5 

ns = l-{nD + JTL + KM + ni), (5) 

where no, til, km, Xi and ng are the proportions of cells that become defectively infected, latently 
infected, persistently infected or activated virus-producing cells or die, respectively. 



2 Alternative models 

The results from the main text were based on our default model but we also constructed 12 alter- 
native models (Table SI) that differ in one or more characteristics of the HIV-1 life cycle: 
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• No CD4 + : We did not include the number of CD4 + T cells per |_d in the model fits. 

• 6i 6 =0.5: The average lifespan of virus-producing cells (J 6 and M 2 ) is longer (l/5/ 6 = 2 
days). 

• e = 1.0 : Drug efficacy is assumed to be 100%. 

• Direct: The fate of an infected cell is determined directly after infection and not sequen- 
tially during the intracellular eclipse phase. 

• No eclipse: There is no intracellular eclipse phase, i.e., no stages h to I5. 

• Horn. 1 : Latently infected cells were assumed to be a homogeneous cell population (no L 2 ) . 

• Horn. 2: Persistently infected cells were assumed to be a homogeneous cell population (no 
Mi). 

• Prod./death 1: Persistently infected cells (M 2 ) and activated, virus-producing cells (I 6 ) have 
a different viral burst size. 

• Prod./death 2: Persistently infected cells have a different death rate {8m 2 = St) than acti- 
vated, virus-producing cells (<5/ 6 ), but the same viral burst size N. 

• Prod./death 3: Persistently infected cells (M 2 ) and activated, virus-producing cells (I 6 ) have 
a both a different viral burst size and death rate than activated, virus-producing cells. 

• Burst 1: Transcriptional bursts in latently and persistently infected cells last for one hour 
on average. 

• Burst 2: Transcriptional bursts in latently and persistently infected cells last for one week 
on average. 

Table SI provides an overview of all models together with the sum of squared residuals (SSR) 
from fitting them to the data of the five patients. The estimated parameters of the HIV- 1 dynamics 
for all models are given in Table S2. 

It can be seen that the default model provides the best fit to the data together with the model 
where the fate of an infected cell is determined directly after infection (Direct). The model that 
does not include the number of CD4 + T cells {No CD4 + ) results in an only slightly smaller SSR 
than the default model. Different assumptions about the death rate of virus-producing cells 
(5/ 6 = 0.5) and the drug efficacy {e = 1.0) do not improve the quality of the fits. Assuming a 
different viral burst size for persistently and activated virus-producing cells {Prod./death 1 and 
Prod./death 1) is also consistent with the data. However, the estimated viral burst sizes in per- 
sistently infected cells were higher than in activated, virus-producing cells. This seems to be a 
rather unrealistic scenario. A longer duration of transcriptional bursts results in a reasonably 
good model fit {Burst 2) whereas bursts with a duration of only one hour are not consistent with 
the data {Burst 1). All other models {No eclipse, Horn. 1, Horn. 2, Prod./death 2 and Prod./death 2) 
fit the data substantially worse than the default model. Note that some models do not converge, 
indicating unrealistic assumptions about the HIV-1 life cycle and/or too many free parameters 
(overfitting). 
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Table SI. Comparison of the alternative models with the default model from the main text. 



Model 


CD4 + 


Cell fate 


Eclipse phase 


Latent 


Persistent 


Production 


Death 


Burst 


No. of parameter 


s Convergence" 


SSR b 


Default 


Yes 


Sequential 


Yes 


Het 


Het 


Same 


Same 


Day 


12 


Yes 


160.7 


No CD4 + 


No 


Sequential 


Yes 


Het 


Het 








11 


Yes 


159.1 


6 h = 0.5 


Yes 


Sequential 


Yes 


Het 


Het 


Same 


Same 


Day 


12 


Yes 


165.6 


e = 1.0 


Yes 


Sequential 


Yes 


Het 


Het 


Same 


Same 


Day 


12 


No 


161.7 


Direct 


Yes 


Direct 


Yes 


Het 


Het 


Same 


Same 


Day 


12 


Yes 


160.2 




Yes 


Direct 


No 




Het 


Same 


Same 


Day 








Horn. 1 


Yes 


Sequential 


Yes 


Horn 


Het 


Same 


Same 


Day 


11 


No 


221.3 


Horn. 2 




Sequential 


Yes 




Horn 


Same 


Same 




11 


No 


274.2 


Prod. /death 1 


Yes 


Sequential 


Yes 


Het 


Het 


Different 


Same 


Day 


13 


No 


157.6 


Prod./death2 


Yes 


Sequential 


Yes 


Het 


Het 


Same 


Different 


Day 


12 


No 


197.6 


Prod./death3 


Yes 


Sequential 


Yes 


Het 


Het 


Different 


Different 


Day 


13 


No 


167.3 



Burst 2 Yes Sequential Yes Het Het Same Same Week 12 Yes 165.5 

a No if the fitting algorithm failed to convergence for one or more patients 
b Sum of squared residuals 



Table S2. Estimated parameters of HIV- 1 dynamics for all models. Values are given as geometric means over all five patients. For dimensions, we refer to Table 
1 from the main text. 



Model 




A 


P 


8j 


S D 


fD 


h 


/m 






a 


JV 




Default 


1.14e-02 


6.84e+01 


1.35e-06 


1.65e-02 


8.25e-05 


1.41e-01 


3.51e-03 


2.50e-01 


7.89e-02 


1.03e-01 


2.67e-03 


2.14e+04 




6 k = 0.5 


8.44e-03 


7.24e+01 


1.07e-06 


2.71e-02 


6.72e-05 


2.28e-01 


7.15e-03 


2.16e-01 


7.27e-02 


1.16e-01 


3.26e-03 


3.16e+04 




e= 1.0 


2.28e-02 


1.01e+02 


4.81e-07 


7.39e-03 


3.00e-04 


1.05e-01 


4.02e-03 


1.74e-01 


6.31e-02 


8.64e-02 


8.37e-03 


2.48e+04 




Direct 


2.06e-02 


6.98e+01 


6.86e-07 


9.32e-03 


1.32e-04 


1.35e-01 


2.74e-03 


1.68e-01 


8.09e-02 


1.70e-01 


2.14e-03 


2.06e+04 




No eclipse 


1.91e-02 


5.46e+01 


6.77e-07 


6.75e-03 


6.99e-05 


1.65e-01 


3.02e-03 


7.51e-02 


6 81c-02 


1.33e-01 


2.40e-03 


2.07e+04 




Horn. 1 


1.61e-02 


6.12e+01 


2.38e-06 


7.79e-03 


5.60e-05 


1.63e-01 


4.41e-03 


9.72e-02 




5.25e-02 


4.94e-04 


7.98e+03 




Horn. 2 6.73e-02 1.19e+02 3.42e-07 3.74e-03 4.14e-04 1.20e-01 1.29e-02 4.87e-03 1.47e-01 3.89e-02 6.08e+03 


Prod. /death 1 


1.99e-02 


1.59e+02 


6.03e-07 


2.63e-02 


4.19e-05 


1.12e-01 


2.77e-03 


5.22e-02 


7.63e-02 


8.62e-02 


5.44e-03 


1.59e+04 


3.73e+04 


Prod./death2 


9.51e-03 


8.54e+01 


2.46e-07 


2.95e-02 


1.44e-04 


2.67e-01 


6.13e-03 


4.96e-02 


1.80e-02 


9.62e-02 


5.73e-04 


1.24e+05 




Prod./death3 


7.94e-03 


1.14e+02 


1.23e-06 


4.73e-02 


1.15e-04 


2.57e-01 


7.53e-03 


7.26e-02 


4.80e-02 


2.99e-02 


2.60e-03 


5.22e+04 


4.49e+03 


Burst 1 


3.89e-02 


9.06e+01 


3.33e-07 


6.40e-03 


3.08e-04 


1.36e-01 


2.65e-03 


3.43e-02 


8.92e-02 


3.08e-01 


1.76e-03 


1.77e+04 




Burst 2 


7.27e-03 


9.98e+01 


3.23e-06 


4.62e-02 


1.38e-03 


1.60e-01 


5.06e-03 


2.17e-01 


1.34e-02 


6.68e-02 


3.80e-03 


1.65e+04 





3 Model fits to data from other patients 

The fits of the default model (see main text) to the data of the four additional patients (103, 104, 
110 and 111) from Fischer et al. [1] are shown in Fig. S1-S4. 
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Figure SI. Model fit for patient 103. Circles denote measured data and lines represent the best model fit. 
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Figure S2. Model fit for patient 104. Circles denote measured data and lines represent the best model fit. 
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Figure S3. Model fit for patient 1 10. Circles denote measured data and lines represent the best model fit. 
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Figure S4. Model fit for patient 111. Circles denote measured data and lines represent the best model fit. 
We did not include the outlier in plasma HIV-RNA at 308 days after start of treatment in the model fit 
because treatment had to be interrupted due to side effects. 
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4 HIV- 1 DNA positive cells 

The study by Schmid et al. [2] measured the number of HIV-1 DNA copies per 10 6 peripheral 
blood mononuclear cells (PBMCs) in patients that initiated cART during the acute {n = 24) and 
chronic infection (n = 15). We assumed that there is only one proviral DNA copy per infected cell 
[3], and scale the number of HIV-1 DNA copies with the average number of PBMCs per ml from 
the five patients in the study by Fischer et al. [1]. This resulted in numbers ofHIV-1 DNA positive 
cells per ml that are expected to exist during acute (mean: 522 cells) and chronic infection (mean: 
1964 cells). Two standard deviations around the mean of those numbers are shown in Figure 4 
from the main text. 
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